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Abstract 

We present a new approach to the static finite temperature correlation functions of 
the Heisenberg chain based on functional equations. An inhomogeneous generalization 
of the n-site density operator is considered. The lattice path integral formulation with a 
finite but arbitrary Trotter number allows to derive a set of discrete functional equations 
with respect to the spectral parameters. We show that these equations yield a unique 
characterisation of the density operator. Our functional equations are a discrete version of 
the reduced q-Knizhnik-Zamolodchikov equations which played a central role in the study 
of the zero temperature case. As a natural result, and independent of the arguments given 
by Jimbo, Miwa, and Smirnov (2009) we prove that the inhomogeneous finite temperature 
correlation functions have the same remarkable structure as for zero temperature: they 
are a sum of products of nearest-neighbor correlators. 



1 Introduction 



The seminal model of integrable lattice systems is the spin- 1/2 Heisenberg chain whose spec- 
trum can be exactly calculated by Bethe ansatz [TJ[2l. Many properties of the model have 
been calculated, e.g. the thermodynamical potential [3j|4j and the scaling dimensions [5|-|11| 
describing the long distance asymptotics of all correlation functions. 

The calculation of space and time dependent correlation functions is a notorious problem. 
A well-known procedure is the form factor approach which is based on a spectral decompo- 
sition of the space and time dependent two-point function by use of the eigenvalues of the 
Hamiltonian and of the matrix elements of the local operators with respect to the eigen- 
states of the Hamiltonian. For interacting systems the calculation of matrix elements is 
nontrivial 12 -14 . In applications the spectral sum is carried out over sectors with restricted 
particle number [15-20 . Recently, analytic summations over infinite classes of states were 
achieved |2lH24[ 



Methodically different are approaches that aim at the computation of the reduced density 

func- 



matrix. Based on for instance the representation theory of quantum algebras 113 25 



tional equations |26|, or the algebraic Bethe ansatz [27-31 , multiple integral representations 
for the elements of the density matrix were derived. 



In a series of works [46j-[50j it was shown that the general ground-state correlation func- 
tions of an inhomogeneous version of the model can be represented by algebraic expressions 
involving only a few independent functions describing nearest-neighbor correlations. 



An exponential formula for the reduced density matrix was obtained in |51[|52| . On 
the basis of an explicit calculation for short distances of local operators, this formula was 



conjectured in 53 -55 to be also valid for finite temperature or finite length. In |56l 57] a 



fermionic structure on the space of local operators was identified from which a generating 



function of all local correlation functions was obtained 58 for a very general inhomogeneous 
vertex model including the finite temperature and the finite length Heisenberg chain as special 
cases. Even in this most general situation the inhomogeneous correlation functions depend 
on only two functions [3lJ[58] . 

The properties of the Heisenberg chain at finite temperature are obtained within the quan- 
tum transfer matrix approach using the correspondence between the Heisenberg chain and 
the six- vertex model. This method has been developed as an alternative approach |32[]43| , 



see also the reviews 44 , 45 , to the thermodynamical Bethe ansatz which is based on combina- 
torics and the controversial string picture of the excited states of the Heisenberg chain. Here 
we consider independent spectral parameters in both the chain and the quantum direction of 
the six- vertex model to define an inhomogeneous n-site density operator as introduced in the 
work [25] . For the zero temperature case, a functional equation of Knizhnik-Zamolodchikov 



type |59[|60| was derived 13,26 . At finite temperature the functional equation is not satisfied 
as it receives finite temperature correction terms which seem unmanageable. 

We find the functional equation approach to the computation of correlation functions very 
important [ 61] . Such an analysis is basically built only on the fundamental properties of the 
i?-matrix: the Yang-Baxter equation, the 'standard initial condition', and also, if satisfied, 
the crossing symmetry. Hence, a generalization of previous work to higher spin or higher 
rank models looks most feasible along this way. Before doing so, it is important to first 
understand the computation of the finite temperature correlations of the Heisenberg chain by 
use of suitable functional equations. 

Facing the above mentioned problem, the first and obvious intention is to identify the finite 
temperature correction terms appearing in the functional equations. In this paper, however, 
we perform an analysis based on a rather contrary reasoning. We find that the correction 
terms vanish at certain discrete points of the spectral parameters. Interestingly, the reduced 
information is still sufficient to fix the full functional dependence of the finite temperature 
correlations. We therefore call these equations 'discrete functional equations'. This paper is 
organized as follows. In section 2 we define the inhomogeneous density matrix and derive 
the discrete functional equations in section 3 with the proof of uniqueness deferred to the 
appendix. The solution to the discrete functional equations for the XXX case is presented 
in detail in section 4. Generalizations to the XX Z case and for systems with an a-seam are 
found in the appendix. 



2 The inhomogeneous n-site density operator 
2.1 The i?-matrix 

The Hamiltonian of the spin- 1/2 Heisenberg chain is given by 



H XXX = £ (of af +1 + a\a y i+l + Aa*a* +1 ) (1) 



i=i 
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where a x , a y , and a z are the Pauli matrices 



a 



1 

1 0/ ' 



a" 



-i 

1 



a" 



1 
-1 



We denote the orthonormal eigenstates of a z by v + (eigenvalue 1) and v~ (eigenvalue -1). 
We are interested in the thermodynamic limit of ([!]). Therefore, the boundary conditions 
are unimportant, but we may impose periodic boundary conditions whenever arguments are 
based on a finite size L. 

The i?-matrix of the classical counterpart of the Heisenberg chain, the six-vertex model, 
is given by 

/o(Ai,A 2 ) \ 
1 6(Ai,A 2 ) 
b(Xi,X 2 ) 1 
V a(Ai,A 2 )/ 

For A = ch(??) the functions a and b are defined as 



£(Ai,A 2 ) := 



with 



a(Ai, A 2 ) = [1 + Ai - A 2 ] g , 6(Ai, A 2 ) = [Ai - \ 2 } q . 
sh(rjx) 



(2) 



sh(r/) 



for q 



In the limit q — > 1 the functions for the XXX case (A = 1) are obtained: 

«(Ai, A 2 ) = 1 + Ai - A 2 , 6(Ai, A 2 ) = Ai - A 2 . 



(3) 



The matrix elements of the i?-matrix are considered as local Boltzmann weights of the six- 
vertex model 

R(\i,\ 2 )(v^ ®v^) = ^(Ai,A 2 )^ 2 2 



0"1,(T2 



and are graphically represented by vertices with spin variables at the ends of the bonds and 
spectral parameters associated with the lines, see Fig. [TJ The i?-matrix enjoys the following 
properties: 



/'i 



/'2 



Atx \ Ai 

(71 <7 2 

Figure 1: Graphical depiction of the matrix element R(Xi, A 2 )ctJct2 • 
Yang-Baxter equation (illustrated in Fig. [2]) 

£i, 2 (A 2 , A 3 ) ^ 2j3 (Ai, A 3 ) i?i, 2 (Ai, A 2 ) = i? 2 , 3 (Ai, A 2 ) i?i, 2 (Ai, A 3 ) ^ 2 , 3 (A 2 , A 3 ), (4) 
initial condition (illustrated in Fig. [3]) 

R(\,X) = I®I, (5) 
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and unitarity (illustrated in Fig. [2J 

R{\ 2 ,\ 1 )R(\ 1 ,\ 2 ) = C V 1(8)/, C u :=a(X 1 -X 2 )a(-(X 1 -X 2 )). (6) 
Furthermore it possesses the crossing symmetry (illustrated in Fig. [3]). For 

4 

i?(A l5 A 2 ) := PR(\ U X 2 ) =J2(<Pi® i>j) 

»J=l 

with suitable maps 6 End(V) and with P denoting the permutation operator on two 

sites we find 

fa ® V>j) = (S ® I) R(X 2 - 1, Ai) (S ® I) , (7) 

»J 

where 5 is defined by 

S(t> + ) = S*^ - ) = — v + . 

A graphical depiction of the matrix elements of S with respect to the basis v + ,v~ is introduced 
in Fig. [4} 



Mi M2 A*2 A*3 Mi M2 




Figure 2: Graphical depiction of the Yang-Baxter equation (left graph) and unitarity (right 
graph). A summation over all allowed values of the spin variables on closed bonds is un- 
derstood implicitly. Ai,A2,A3 denote spectral parameters assigned to lines; cri,a 2 ,as and 
Mi)M2)M3 are s Pi n variables assigned to open bonds. 




Figure 3: Graphical depiction of initial condition (left graph) and crossing symmetry (right 
graph) . 
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Figure 4: Matrix elements of the identity and the single site map S. 



2.2 Construction of the inhomogeneous n-site density operator 

We study a semi-infinite six-vertex model on a regular lattice of height N and infinite width, 
see Fig. [5j The width corresponds to the length of the Heisenberg chain, the height corre- 
sponds to the Trotter number of the imaginary time discretization. For a suitable choice of 
spectral parameters assigned to the rows, the product of the N row-to-row transfer matrices 
approaches exp(— f3H) in the limit N — > oo. The following analysis, however, is done for finite 
N and the limit N — > oo is taken at the very last stage. 

A very important object for the study of the model on the semi-infinite cylinder is the 
column-to-column transfer matrix. This matrix (in the limit ./V — > oo) is also known as the 
quantum transfer matrix of the Heisenberg chain. We define the corresponding inhomogeneous 
quantum monodromy matrix as 

T(A) := R N)N+1 (X, u N )--- # 2 ,3(A, f2)-Ri,2(A, 

With respect to the basis v + ,v~ in the space V\ the matrix is denoted as: 



r+(A) r+(A) 

T+(A) 7T(A) 



+m r-m). 7^(A) G End (Vq) , V Q := V U1 ® • • • ® V m 



The quantum transfer matrix is obtained by taking the trace over the auxiliary space Vy. 

T(A) := try A (T(A)) = T+(A) + 7T(A) 
We define the inhomogeneous n-site density operator 

D n {X u . . . , A n ) G End(Vxi ® • • • ® VxJ 

in the following way: 

d fA a v- i™ tr VQ([ r (°)] fc nAi)r(A 2 )--r(A w )[r(o)]*) 

' ' • ' nj fe™ tr yQ ([T(0)]^r(A 1 )T(A 2 ) • • • T(A„)[T(0)]*) l8j 

Due to the commutativity of the quantum transfer matrices, [T(A),T(//)] = 0, the eigenvec- 
tors of T(A) can be chosen independently of the spectral parameter A. We assume that the 
spectral parameters v\, . . . , i/jy are chosen in such a way, that the leading eigenvalue of T(0) 
is non-degenerate. In this case only the (normalised) eigenvector &o, corresponding to the 
leading eigenvalue Ao, contributes to the trace: 

n x ^ ($o|T(A 1 )r(A 2 )---r(A n )|«l>o) 

lli=l A o(Aj) 



Note that for the physical values stated in equation (11) there is indeed a finite gap between 



the leading eigenvalue and the rest of the spectrum even in the Trotter limit iV — > oo. For a 
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Figure 5: Graphical depiction of tr Vq {[T{0)} k Ta^ {X^T^ 2 {X 2 ) • • • T£ l (X n )[T(0)] k ) 



proof see 35,36 and chapter 13 of the book 62 where an intuitive argument based on the 



high temperature case is given. 

The concatenation of transfer respectively monodromy matrices in formulae ^ , ^ is to 
be understood with respect to the quantum space V Ul <g> • • • <S> V VN , meaning that the matrix 
elements are given by 



($o| V l 1 1 (X l )T^(X 2 ) • • • 7£"(A W ) |$p 
nr=iAo(A,) 



(10) 



The re-site density operator for the infinite Heisenberg chain in a heat bath with temperature 
T is obtained from D n (Xi, . . . , A n ) by first specialising the parameters < j < iV in the 
quantum space (for even Trotter number) to the values 



P 



8 , N 

•£ + 1 * = l,...y, 



(11) 



then performing the Trotter limit A r — )• oo and taking the so-called homogeneous limit Aj — > 
for 1 < i < n. 



3 The discrete functional equations 

For zero temperature the density operator satisfies the functional equation 
A n (Xi, . . . , A n ) [D n (Xi, . . . , A n _i, A n )j = D n {X\, . . . , A n _i, X n - 



1) 



fl2al 



for arbitrary complex A n . Under suitable asymptotic conditions and analyticity the solution 
to this functional equation is unique |49|50|. At finite temperature resp. finite Trotter number 



the functional equation does not hold for arbitrary A n . The proof of (12a) for T = is based 
on so-called 'Z-invariance' arguments [2] relating the partition functions of different lattices. 
At T = certain terms disappear as they can be viewed as boundary terms at infinity. For 
finite Trotter number these terms matter. When we attempted to quantify these terms we 
were surprised to realize that the unwanted terms do not appear for the finite set of arguments 



A n = V\, . . . , UN. 



(12b) 
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Interestingly, equation (12a) for (12b) and asymptotic and analytic properties of the density 
operator are sufficient to uniquely fix the solution 161] as will be shown at the end of this 
section. 




Figure 6: Graphical illustration of the matrix element (A n [D n (\i, . . . , A n )])^"'^™ which is 
basically given by the partition function of a lattice as shown in Fig. [5] extended by a sequence 
of i?-matrices. Here, the infinitely many columns in the thermodynamical limit have been 
replaced by the leading left and right eigenstates <3?o- Note that for the normalization of the 
density matrix we would have to multiply it with C^ 1 Yli=i Ao(Aj). 




Figure 7: Graphical illustration of (A n [D n+1 (Xi, . . . , A„, A„ - l)])^'."^^ 1 multiplied by 
C^ 1 [Y\i = i Ao(Aj)] Ao(A n — 1). Note that a contraction of the 'A n -line' will reduce the ob- 
ject to a n-site density matrix. The same holds true for a contraction of the 'A n — 1-line'. 

The operator 

A n (X u . . . , A n ) e End (End(y Al ® • • • ® V X J) 
is constructed from i2-matrices, see Fig.pl The image of an operator B 6 End(V® n ) is defined 



7 



a£] 

A„(Al, . . . , X n )[B] := Ca tTv„ [Rl,2(M, An) • • • -Rn-l,n(A n -l, A„) (Pr s )n,n+1 

{B ® I n+ \) i? n _i jn (A n , A n _l) ••• i?2,3(An, A2) -Rl,2(A n , Ai 

The constant is defined as 

n-l 

Ca := C^ X (A„, A;) 



(13) 



and Pr, 



s } ( s I denotes the (unnormalised) projector onto the 2-site s^-singlet 

s = v + (g) f ~ — v ~ f + . 



(14) 



In order to derive the discrete functional equations (12a) for (12b) we consider the inhomo- 
geneous density operator on (n + 1) sites with the specialization A n+ i = X n 



-Dn+i(Ai, . . . , A n , A n — 1) 



($o| T"(Ai)T(A 2 ) • • • T(\n)T(X n 



1: 

l)l*o> 



Ao(Ai)A (A 2 )---Ao(A n )A (A n -l) 



We consider the action of A n on the first n sites of D n+ i(\i, . . . , A n , X n — 1), see figure Q. 
Taking the trace of the image ^4„[-D n+ i(Ai, . . . , A n , A n — 1)] with respect to the n-th respecively 
(n + l)-th space yields 



try n (A n [L» n+1 (Ai, . . . , A n , A n - 1)]) = D n (Xi, . . . , A n - 1), 
try n+1 (A n [D n+ i(Xi, . . . , A„, A„ - 1)]) = A n [D n (Xi, ... , X n )] 



(15) 
(16) 



for all values of X n . For the special case that A ra is equal to one of the horizontal inhomoge- 
nieties fx, . . . , vp?, an important simplification occurs due to the standard initial condition 
(equation (IS]) ) leading to the factorization of two i?-matrices. Consequently the left hand 



sides of equations (15) and (16) become equal and hence the right hand sides are identical. 
In other words we obtain equation (12a). The graphical proof is given in Fig. [8j [^] 

Note that the presented construction can be substantially generalized to other boundary 
conditions. If (<3?o| an d |^o) m ([£]) an d Fig. [6] etc. are replaced by arbitrary left and right 
eigenstates ($/| and |<F r ) (to the same or different eigenvalues of the quantum transfer matrix) 



all arguments leading to (15), (16) remain valid cum grano salis, see also [63]. 



3.1 Characterization of D n 

Note that the functional equations can be formulated for every spectral parameter Aj using 
the symmetry 

Ri,i+i(X-i, Xi+\)D n (. . . , Aj, Aj+i, . . . )iij j i-|-i(Aj+i, A») = D n (. . . , Aj+i, Aj, . . . ). 

However, to give a (recursive) characterization of D n (Xi, . . . , X n ) it suffices to regard the 
functional equations for the last spectral parameter alone: 



The definition of A n is more eas ily st ated by using the duality End(V) — V ® V* . (Compare eq. (341.) 



It is possible to state equation (12a I for a more general model, including a magnetic field and disorder 



parameter. This issue is discussed in appendix [C| 
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Figure 8: LHS: Using the crossing symmetry for the quantum monodromy matrix T{X n — 1) 
acting on site n + 1 the direction of the (A n — l)-line is changed, thereby raising the spectral 
parameter by one. For the special case that A n is equal to one of the parameters Vi (we choose 
v\) we find the RHS by use of the initial condition and unitarity. RHS: the two traces in 



eqs. (15) and (16), contractions over the pair (a n ,fi n ) and (a n +i, Mn+i) ; respectively, yield 



the same object. 



The n-site density operator D n (Xi, . . . , A n ) is uniquely characterized by the properties: 

A n (Xi, . . .,Xn-i,Ui)[D n (Xi, . . . , A n _i, v^] = D n (X 1 , . . .,X n -i,i>i - 1) i = 1...,N (17) 
and 

lim D n (Ai,...,A n ) = -D n _i(Ai,...,A„_i) ®I (18) 

for generic values of the inhomogeneities i>i. The proof of this theorem is based on the 
functional dependence of the operator L> n (Ai, . . . , A n ) on the n-th spectral parameter. For 



the XXX case each matrix element ( 10 ) is a n-variate polynomial in the variables Aj of degree 
at most N divided by a (known) polynomial of the same type. For the XXZ case the same 
statement holds with respect to polynonials in the variables q 2Xi . As a consequence it suffices 



to show that the N + 1 conditions given by equations (17) and (18) are independent. We 
state the proof in appendix |B} 

4 The n-site density operator for the XXX chain 

For the XXX case we discuss the explicit solution for the inhomogeneous n-site density oper- 
ator for finite Trotter number in full detail. The solution takes the same form as in the zero 
temperature case where only the function oj has to be replaced by its finite Trotter number 
counterpart. More precisely, we will show that the solution for arbitrary n can be stated in 
terms of a single nearest neighbor correlator oj in the form 

[n/2] / rn \ 

D n (X 1 ,X 2 ,...,X n ) = J2J2\ INV^) /n,/,j(Ai,A 2 ,...,A„). (19) 

m=0 I,J \p=l I 



9 



The structure coefficients f n l J £ End (V® n ) are the same as constructed in 49 for the zero 



temperature case, quite explicit formulas can be found in 164 65 . These coefficients are some 



times called the 'algebraic part' of the density matrix. The nearest neighbor correlator oj is 
called the 'physical part' which aquires a temperature dependence. This simple correspon- 



dence of zero and finite temperature was conjectured in 53 where the explicit factorization 



of multiple integral formulas for short distance correlators was carried out. Here we complete 



our proof of (19) on the basis of discrete functional equations which is alternative to the 



algebraic proof by Jimbo, Miwa, and Smirnov |58|. 



4.1 The 2-site function oj 

In the single site case the density matrix is completely fixed by symmetry, in the two site case 
we find a linear combination of the identity operator and the projector Pr s onto the singlet 



state (14) 



Di(Ai) = iid, £> 2 (Ai, A 2 ) = i id® id + | w(Ai, A 2 ) Q id® id- Pr s ) . (20) 



Equation (20) yields 



tr (D 2 (m, A) Pr s ) = ~-u(n,\). 



(21) 



From the properties of Z?2(Ai, A 2 ) we obtain important properties of oj which are needed for 



the proof of equation ( 19 ) 



1) The function oj obeys the discrete functional equation^] 

3 



(\ \ -n ^2i(A 2 i - 2) 
w(Ai, A 2 - 1) = — -2 w(Ai, A 2 ) 



1 - A 



21 



2 (1 - A10 



for A 2 = vi,...,v N . (22) 



2) The asymptotic behaviour is given by 



^(Ai, A 2 ) > 0. 

A2— >oo 



(23) 



3) oj is symmetric under exchange of arguments: 

w(Ai,A 2 ) = w(A 2 , Ai). 



(24) 



4) Furthermore Ao(Ai)Ao(A 2 ) w(Ai, A 2 ) is a 2 variate polynomial of degree iV with respect 
to Ai, A 2 . 

An explicit expression of the function oj in terms of integral equations has been obtained 
in [53] . Here we present a derivation having the advantage that all occurring functions can 
be obtained in terms of a single auxiliary function. To this end, we introduce the modified 
quantum transfer matrix in a larger Hilbert space 



Vq := V Ul (8) K 



;/ 2 



V UN (8) V^+i+e (8) Vp. 



3 We use the abbreviation Aij := — Aj 
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Explicitly 



T(A) := try A (R N+2 ,N+3(^, fi) -RiV+l,Ar + 2(A, /x + 1 + e) R N;N+1 (X, V N ) ■■■ -Rl, 2 (A, Ui)) . 

Compared to T(A) the Trotter number of T(A) is enhanced by two. The parameter e is 
considered as a small deformation. For e = the leading eigenvector $o of T(A) corresponds 
to the tensor product of the leading eigenvector <I>o of T(A) in the space V Vl ® • • • ® K/jv with 



the 2-site s/2-singlet (14) in the two additional quantum spaces. The logarithmic derivative 
of the leading eigenvalue Ao(A) of T(A) with respect to the small deformation parameter e 
evaluates to 



de 



ln(A (A)) 



e=0 



(*0 


( 






(*0 


f(X) 





e=0 



lim 

k— ¥00 



tXy Q lT(ji)[T(p) [£T(\)) T(0) 



e=0 



tr^r(//)|T(0)j T(X)[T(0) 

^v,®v x [D 2 (p, A) R(X, fi) (£R{fi + e, A)) | £=0 
tT v„®v x {D 2 {fi, A) R(X, fi) R(fi, A)) 



(25) 



Equation (25) connects the logarithmic derivative of the leading eigenvalue Ao(A) to the ex- 
pectation value of the operator R(X, fi) [4^R{fi + £, A)) | e= o with respect to the (unmodified) 
density operator D 2 (fi, A). A graphical illustration of equation (25) is given in Fig. [9j 
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Figure 9: Graphical illustration of the last equation in (25), see also Fig. 10 



From definition (21 ) and tr (I^O-f, A)) = 1 we find 
d 



de 



ln(A (A)) | e=0 



1 



, _ (A _^)2 i(1 + A-/x)--+o;(m,A)). 



(26) 
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Figure 10: The i?-matrix symbol with central dot denotes the derivative (^-R(m + A)) | £ =o- 



For the leading eigenvalue Ao(A) the integral expression 

ln(A (A)) = ln(A - \i + 1) + ln(A - \i - e - 1) 

[N/2] N 
+ ^ In (A - u k + 1) + J] ln(A - u k ) 

k=l fc=[jV/2]+l 

-— / 1 

2mJ L (X - - uj + 1) 



ln(l + a(u))du (27) 



can be derived from the Bethe ansatz equations |40[|41| . The auxiliary function a is determined 
by the integral equation 

M i( ^ lo (^i) + ,(^ 



k=i 



fc=[JV/2]+l 
1 



A - v k 



+ 



2vri 



((A-w)3-l) 



In (1 + a(w)) dw. (28) 



The contour L encloses the Bethe roots on the imaginary axis and also the parameters \i and 
Ui with i < [N/2], which we assume to be real and close to zero. [Of course, by multiplying 
the spectral parameters by the imaginary unit i the expressions can be written in terms of 



integrals along the real axis.] Setting e = in equations (27) and (28) yields the description 
for the leading eigenvalue of the unmodified quantum transfer matrix, 

[N/2] N 

ln(A (A)) = ln((A-M) 2 -l)+ £ In (A - i/ fc + 1) + £ m(A - v k ) 

k=[N/2]+l 

1 ln(l + a{u)))duj, 



k=l 



2m J L (A - uj)(\ -u + 1) 
with auxiliary function a(A) = a(A)i determined by 



WW / \ i\ N / \ : 



fc=i 



fe=[JV/2]+l 



\-v k 



+ 



2vri 7 L ((A- W )2-l) 



In (1 + a(w)) dw. (29) 
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From equations (26) and ( |27[ ) we obtain 

1 (1-(A-/,) 2 ) 



duj. 



(30) 



2 2vri J L (A-w)(A-w-l)l + o(w) 

The function G introduced in [29] is here obtained in terms of the modified auxiliary function 



G(A,/i) :=-f ln(a(A))| e=0 . 
ae 



The linear integral equation 
G(A,/x) 



■1 



1 



d<jj 



(31) 



(A - /x)(A - /x + 1) 2vri J L ((A - cj) 2 - 1) 1 + a(u;) 
follows from equation (|28|). Independent of our construction it is possible to take equations 



(29)-(31) as definitions and read off the desired polynomial properties as well as (22)-(24). 



4.2 The algebraic part 



We adopt the notations of 49,50 to describe shortly the construction of the algebraic part. 

(32) 



It is convenient to consider the vector 

/i n (Al, • • • , A„) := (D n (Xi, . . . , A„) ® 7® n ) (s n ), 
with the following numbering of the spaces: 

hn(Xi, . . . , A n ) G Vi <8 V 2 <8 • • • <8 F„ <8 Vk (8) • • • <8 ® Vj 
The vector s n G y® 2n i s defined by nested 2-site sfo singlets 



n 

i=i 



Sj j, with s = v 



v — v <S> v 



(33) 



Note that under the duality (32) the vector s n corresponds to the identity operator on n 



sites. Using the duality transformation and the crossing symmetry, A n E End (y® 2n ) takes 
the form 

4° = (]I ^(Al, AO J (-l) n i?l,2(Al,2 " 1) • • • Rl,n(\l,n ~ l)i\ii?l,n(Al,n) ■ ■ ■ Rlfl(\l?). 

(34) 

In this section we consider the functional equation with respect to Ai. In |49| a family of 
operators 



nX^'(\l, . . • , An) € Horn ( V ^n-2)^2n\ j 1 < i<J < 



n, 



called X-operators in the sequel, is constructed via transfer matrices whose dimension of the 
auxiliary space has been continuously extended to C [A] . By construction the action of A n on 
the X-operators is given by 



4 1} (Ai,...,A„) 

(Al, ■ • • , A n ) 
f —nX^l(\\ — 1, A2, • • • , A n ) 



for i = l, 



Xn-2&i - 1) A 2 , • • • , A n ) A^_ 2 (Xi, . . . , Ai, . . . , Xj, . . . , A n ) else. 



(35) 



13 



The symbol Aj denotes the missing of the spectral parameter Aj in the sequence of arguments. 
The image of s n evaluates to 



4 1 )(A l! ...,A n )(s n ) 

= s n - 6A«(A 1 , . . . , A n ) (j^ (C 1 " A fi) A ^i - 2)) _1 n4-2(Ai, • • • , A„)(s n _ 2 ) 
Furthermore, for pairwise distinct indices i,j,k,l the commutation relation 



(36) 



(hi), 



n - * „_2 (Al, • • • i An) n-2-^„_4 (Al, . . . , Aj 



\, ■ . . . , Aj, . . . , A n ) 

1 n.-2^n_4 (Ai, . . . , Afc, . . . , A/, . . . , A n ), (37) 

holds, where k' denotes the position of A& in the sequence Ai, . . . , Aj, . . . , Xj, . . . , A n . 



n^n^2 (Alj ■ • • ) A n ) 



4.3 Solution for the n-site density operator 



The X-operators and the 2-site correlation function w, discussed in section 4.1 are joined in 
the definition of fi-operators: 



n {i ' j) - 



1 



~2 rw(Aj, Xj) n X^l(Xi, . . . , A n ) 



The result for h n is given by 



n/2] 



h, 



( (Al, . . . , An) - ^2 L_2m 







A-l,(ilJl) ^K 2 ,(i 2 ,j 2 ) ••• ^A"m,(im,im)( S "-2m)- (38) 



m=0 7,J 

For X = {ki, . . . , /c m } with 1 < ki < ■ ■ ■ < k m < n the shorthand notation 



&K,{ki,kj) : —m ^m'-2(Afci> • ■ ■ ) Afc m ) 



has been used. The inner sum on the r.h.s. of (38) is performed over pairs of m-tuples, I = 
. . . ,i m ), J = (ji, . . . , j m ) of distinct elements of K\ = {1, 2, . . . , n} under the condition 
%i < ■ ■ ■ < i m and i\ < ji, . . . , i m < j m . The subsets are defined as K p+ \ := K p \ {i p ,j p } for 
p > 1. This expression for the inhomogeneous n-site density operator takes the same form as 
for the case of zero temperature. The only change is the replacement of the nearest neighbor 
correlator by the analogue function for the case of finite Trotter number. (This function u 
determines the analytic properties of the expression.) 

In order to prove that (38) provides indeed the re-site density operator it is sufficient to 
check whether this expression is subject to the characterisation derived in section 3.1 From 
equations (35)-(37) for the X-operators and the discrete functional equations (22) for u) it 
is straightforward to verify that the discrete functional equations (17) for the n-site density 



operator are fulfilled. To prove the asymptotic condition (18) one further needs 49 
1 



lim 



nX^liXt, . . . , An)(s„_2) = const. 



A, (X'f: - 1) 



X^(Al,...,An)( S „- 2 ) 



» Si T „-2^n-4(A2, • • • , A n )(s„_ 4 ),i / 1 



Ai— >oo 
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The essential part of the proof that the r.h.s. of (38) multiplied by YYi=i Ao(Aj) is indeed an 
n-variate polynomial of degree N is to show that the additional poles coming along with the 
A-operators vanish. This is achieved along the same lines as in |50| by rewriting equation (38 ) 
in exponential form and using the symmetry property of u; ( 24 ) to show pairwise cancellation 
of residua. 



5 Discussion 

In this paper we showed how to derive the correlation functions of six-vertex models on 
a lattice with infinite width and arbitrary height N. These systems comprise as a special 
case, in the limit of infinite Trotter number N — > oo, the spin-1/2 Heisenberg chain in the 
thermodynamic limit at finite temperature. Also, but not explicitly treated here, the general 
six- vertex model comprises the zero temperature Heisenberg chain on a finite ring |54||61| . 

The analysis was based on discrete functional equations (of Knizhnik-Zamolodchikov 
type). For a certain set of N many values of the spectral parameter we showed the va- 
lidity of the functional equations which for zero temperature and field are known to hold for 
all complex values. Despite the restricted validity, the discrete functional equations suffice to 
uniquely identify the solution. For finite Trotter number N, the dependence on the spectral 
parameter is polynomial of degree N. A formal proof of the uniqueness was given. 

A similar treatment of the correlation functions is possible for the anisotropic spin-1/2 
Heisenberg chain with finite temperature as well as magnetic field, and even with a magnetic 
seam. Applications to related Temperley-Lieb models and integrable higher spin-S Heisenberg 
chains |66, 67] are currently under investigation. 
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A Solution for the n-site density operator of the XXZ chain 

For the XXZ case two independent nearest neighbor correlators u) and ui are necessary to 
state the result for D n (X\, . . . , A n ). These two functions can be defined as follows 

w(Ai, A 2 ) := 1 (\ ch(r?) - (s(q)\ D 2 (X l ,X 2 ) \s(q))) , 



w(Ai, A 2 ) := (<a(g)| D 2 (X 1 , A 2 ) \s(q)*) - < S (g)*| D 2 (A 1; A 2 ) \s(q)}) 



ch(r/Ai 2 ) sh(?7Ai 2 ) 

W(Ai, A 2 J 



ch(ry) sh(r/) 
Here we used the notation 

s(q) '■= q~ 1/2 v + ® v~ - q l/2 v~ ® v + , s(qf := q~ 1/2 v~ ® v + - q l/2 v + ® v~. 
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The functions are the finite Trotter number analogues of the functions used in [50] . They are 
subject to the discrete functional equations 

w(Ai,A 2 - l)+w(Ai,A 2 ) +p(Ai,A 2 ) = 

w(Ai, A 2 - 1) + w(Ai, A 2 ) + o;(Ai, A 2 ) + p(Ai, A 2 ) + p(Ai, A 2 ) = 



at A 2 = v\ , 



p(Ai, A 2 ) 
p(Ai,A 2 ) 



1 



[2], 



2 [A12 - l] 9 [Ai 2 - 2] q 4 [A12 - 1]JAi 2 + l} q ' 

3 1 1 [3] 9 

"4 [Ai 2 ] 9 [Ai 2 - 1], 4 [Ai 2 - 2],[Ai2 + 1], ' 



The functions u and Co are symmetric respectively antisymmetric with respect to exchange of 
arguments: 

w(Ai, A 2 ) = w(A 2 , Ai), w(A ls A 2 ) = -u)(A 2 , Ai). 



The vanishing of asymptotics follows from 

(s(q)\D 2 (X 1 ,\ 2 )\s(q)) 

respectively 



1 



A2— >oo 2 



> - ch(r/) 



<«(9)|D2(Ai,A 2 )|s(9)*) 



and (s(g) t | J D 2 (Ai,A 2 )|s('7)) ► -1- 

A2— >oo A2— >oo 



The X-operators constructed in 50 take the form 

n^ n ^ 2 (Al, • • • , A n ) = — Ajj n^n-2(Cl) • ■ • > Cn) +n G n f_ 2 {Clj • • • i Cn)> 



with rational functions n G^\ and n G^2 with respect to the variables Ci := q Xi . The in- 



homogeneous n-site density operator for the XXZ chain is given by (38) with the following 
definition of the ^-operators: 



n^n-2 '•— A?) n^n- 2 (Cl; • • • > Cn) + Aj) n ^i- 2 (Cl> ■ • • ; Cn)- 



B Uniqueness of polynomial solutions to (17) and (18) 



This appendix provides the proof that equations ( 17 ) and ( 18 ) uniquely determine the operator 



D n (Xi, . . . , A n ) for generic values of the inhomogeneities v\,...,vn under the condition that 
the solution for n — 1 sites is known. Recall from section |2.2| that the operator 



Ao(A n )D n (Ai, . . . , A n ) 



(39) 



is a polynomial of degree N with respect to the argument A n for the XXX chain. The 
coefficients of this polynomial are linear operators on V® n , whose matrix entries are rational 
functions of Ai, . . . , A n _i. For the XXZ chain the dependence of p9| ) on A n is of the form 
•^(An) with 

F N (x) := {q~ Nx P(q 2x ) | P polynomial of degree N} . (40) 
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We present the proof of uniqueness for the XXZ case. For pairwise distinct i>{ a basis of the 
space Tn(x) is given by (po>Pi, ■ ■ ■ ,Pn) defined by 



N 



N 



PO 



{*)•■= us? - "t* > i'i { - r): " r s " n 

i=l k=l,k£j '■ 



x - v k \ 



for j = l,...,JV. 



Let Di be the coefficients of the operator (39) with respect to the basis (po, ■ ■ ■ ,Pn) and Di 
the coefficients with respect to a second basis (po, . . . ,Pn) defined by 

Pj(x) := pj(x + 1) for j = 0, . . . , N. 

By construction we have 

Di = ka(vi)D{vi) and A = A (^ - l)D{vi - 1) 

for % = 1, . . . , N. The idea is to regard the coefficients Di as unknown quantities, to obtain a 
system of linear equations for them, and finally to show that this system has a unique solution 
in the generic case. The functional equations (17) are equivalent to the equations 

AoOi 



A n {ui) [Di] 



Ao(i/i - 1 



-Di for i = l,...,N. 



The asymptotic condition (18) fixes the coefficients Dq and Dq: 

D = D = D n -i(\i, . . . , A„_i) ® id. 



(41) 



(42) 



The subspace is spanned by (p\, . . . ,pn) respectively (pi, . . . ,pn)- The elements of the 

matrix Wjv defined by 



N 



Pi = ^( W N)jiPj, 
3=0 



arc 



N 



(W N ) 0J = q~ N 5 0il for l = Q,...,N, (W N ) kfi = ]J[u k - Vj - l} q for k = l,...,N, 

3=1 



and 



(W N ) kl = TT Wk V A ^ for k,l = l,...N. 

3=iU [VI ~ V ^ 



Equations (41) yield a system of linear equations for the coefficients D^. 

N 



A (t/j - 1) 
A (i/j) 



A n {yi) 



3=0 



Di for i = l,...,N. 



(43) 



We choose an arbitrary basis in the vector space End(y® n ) and denote by Dj the coordinate 
vector of Di and by A(z/j) the describing matrix of the linear map ^4 n (fj). Equations (43) 
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yield a system of linear equations for the vectors Dj, i = 1, . . . , N, where the constant term is 
determined by the (known) vector Do- The corresponding homogeneous linear system reads 

D = LAWB. (44) 
The N x N block matrices L , A, and W with blocks of size n 2 x n 2 have the entries 

(W) M :=(^) M 1, (L) M .:=M|=il^l, 

and 

(A) id := 5 id A(ui) for i, j = 1, . . . ,N. 
The solution of ( 44 ) is unique iff 

det (LAW - 1) / 0. 

We make the ansatz 

vx-i =ni + d+-, u 2 i=ni-d+^ for 1 = 1,..., [N/2] (45) 
with a constant parameter d. The Bethe roots n\ are determined by the Bethe ansatz equations 



Jt [n, - v, + 1]„ ™ In, - n k + 1L 



for the leading eigenvalue 



JL IWl, A + 1] iv [«/2 1[A _ , 

3=1 i=l L J J 9 3=1 Z=l 1 JJ <? 

We regard the limit 

n; — > oo under the condition |n; — ^fel^ — )• oo V l,k. 

In this limit the blocks on the diagonal of A are constant and have the only eigenvalues ±1. 
Furthermore WL commutes with A. It thus suffices to show that all eigenvalues of WL are 
different from ±1. We find 

Aq(^,_i-1) [1 + H Ao(na-l) [1-H 



Ao(^z-i) [1 — 2d]_ ' A (^) [1 + 24 



>q 
and 

( ^, 2i _ 1 , 21 _ 1 = !^k 8 (^ )a _ 1 , a = I ^-<,-' M - 1 » 

for /* = 1, . . . , [N/2]. The remaining elements of Wn vanish. The parameter d can easily be 
chosen in such a way that the eigenvalues of the 2x2 matrices on the diagonal of W L are 
not ±1. 
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C Functional equations in the presence of a magnetic field and 
disorder 



The functional equations (12a) for the inhomogeneous n-site density operator of the Heisen- 



berg chain ([!]) hold also in the presence of a magnetic field h, which for the 6- vertex model 

corresponds to a horizontal seam (eW^f FortheXXZ case (A / 1) the introduction of 
a disorder parameter a brings forth a fermionic structure on the space of local operators 56|57| 
and consequently a structural simplification of the solution. For our approach, the additional 
interaction terms do not present fundamental obstacles. Discrete functional equations can be 
derived and analytical properties established. The solution to these equations is again unique. 



Ml 



Ai 




q 2a 
q~ 2a 



,0h , 
e 



Figure 11: Graphical illustration of the matrix element (A n [D n (Xi, . . . , A n )])^"'^™. For the 

normalization of the density matrix we would have to multiply it with YYi=i ^-0* (AO an( ^ the 
inverse of the constant (46). 



The generalized inhomogeneous n-site density operator in the presence of a magnetic field 
and disorder parameter is defined as 



D n (X\, . . . , A r , 



: lim 

k— ¥oo 



try Q ([r( 2 ^+^)(o)]V( /l «(A 1 )r^)(A 2 )---r^)(A n ) [T(^)(o)] fe 

try Q ([T( 2 ^+^)(0)] fc T(^)(Ai)T(^)(A 2 )---T(^)(A„) [T(^)(0)]' 



(2arj+hp) 




T (^)(A 1 )T^)(A 2 )---r^)(A n ) 



/ 



A^ ) (A 1 )A ^)(A 2 )..-Ao^)(A„) 



where 3>q^ denotes the eigenvector corresponding to the leading eigenvalue Aq of the twisted 
transfer matrix 



T (K) (A) :=try A (r N}N+1 (X,v n ) ■ ■ ■ R lj2 {\, u^' 2 ^ 



K G 
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The discrete functional equations hold for the generalized re-site density operator, if the 
definition of the operator A n includes the disorder field in the way shown in figure (11) and 
the constant Ca is set to 



C A :-- 



A, 



AT(An) 



n-1 



(46) 



H) \ A n) 1=1 

Here we do not present the solution to the functional equations but like to refer the reader 
to |3ll[58 . 
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